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ON NEUTRON STAR STRUCTURE AND THE MILLISECOND PULSAR 

% 

Alice K. Harding 

Laboratory for Hi qh Energy Astrophysics 
NASA Goddard Space Flight Center 

ABSTRACT 

The recently discovered millisecond pulsar* (PSR1937+214) is observed to 
be rotating close to the limit of dynamical instability for a neutron star. 

In spite of its extremely rapid rotation, the present measurements of the 
period derivative 2 * 3 put a stringent upper limit on the energy loss from 
gravi tational radiation, thus requiring that the quadrupole moment be quite 
small. The pulsar must also be rotating below the critical frequency at which 
its equilibrium configuration would become non-axi symmetric, since the 
lifetime of this configuration against decay by gravitational radiation is 
very short. This critical frequency, given by the theory of rotating 
ellipsoids^, imposes a restriction on the rotation rate more severe than the 
break-up frequency and may be used to set a lower limit, <p> > 2 X 10*^ g cirf 3 , 
on the density of the star. If the mass is 0.5-1. 5 M 0 , several of the stiffer 
neutron star equations of state may be ruled out, and the radius should be 
less than 16 km. The condition for axi symmetry also imposes an upper limit on 
the rotation rate to which neutron stars may be spun up by accretion disks in 
binary systems, a model recently proposed for the evolution of the millisecond 
pul sar 5 . 


To be published in Nature. 
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Uniformly rotating masses deform into non-spherical equilibrium shapes, 
whose geometry becomes increasingly complex at rapid rotation rates. A 
homogeneous, uniformly rotating body will deform ( in the Newtonian 
approximation) into an oblate spheroid with eccentricity dependent only on 
rotation rate n and density p. At high rotation rates, the equilibrium figure 
can become tri axial at the point where the Jacobi ellipsoids, having lower n 
at a given eccentricity, bifurcate from the axi symmetric Maclaurin sequence. 
This bifurcation point, for a homogeneous body, is^ 


Vil = 0.374 

TTbp 


( 1 ) 


Above this point, and below the bifurcation points of more complex ellipsoids, 
there are actually two possible equilibrium configurations with different n at 
a given eccentricity. It has been shown, however, that the Maclaurin 
spheroids are secularly unstable from either viscosity or gravitational 
radiation reaction at the bifurcation point. 6 Therefore, masses rotating with 
n > n cn * t are expected to become Jacobi ellipsoids. Because of its quadrupole 
moment, though, a Jacobi ellipsoid will relax by gravitational radiation to a 
Maclaurin spheroid with n = ^ C rit ^ spinning up (but losing angular momentum 
by decreasing its moment of inertia) on a time scale of ~0.1 sJ 

There is reason to believe that rapidly rotating neutron stars may deform 
into tri axial ellipsoids. Their rigidity is high and, with the exception of 
possible (small) slippage between the superfluid core and the crust, the star 
is thought to be uniformly rotating. Homogeneity is probably a good 
approximation for neutron stars with masses qreater than around 0.5 M 0 , whose 
density is fairly uniform from the center out to the last few meters of 
crust. Low mass neutron stars have large radii, thick crusts and density 
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profiles which are better modeled by polytropes of Index n > 1. Equilibrium 
figures for poly tropic models have been investigated 8 , and It Is shown that 
Jacobi ellipsoids exist for a poly tropic Index n < 0.8, which would correspond 
to neutron stars of higher mass. By assigning polytropic indices to neutron 
star models, Tsuruta and Cameron 8 conclude that neutron stars with 
M > 0.5 M g can deform into triaxlal ellipsoids. If, at some staqe in its 
evolution, a pulsar is spinning rapidly enough to become non-axi symmetric, 
however, it would not have much chance of being observed, because of the short 
lifetime of the Jacobi ellipsoid phase. Therefore, may be considered an 

upper limit to observable pulsar rotation rates, unless pulsars are very low 
mass neutron stars. 

The millisecond pulsar, with a period of P = 1.558 ms. (n * 4033 s" 1 ), is 
the first observed example of a neutron star spinninq rapidly enough to 
approach the Jacobi bifurcation point and thus affords the possibility of 
constraining neutron star physics. From the presently measured value of the 
period derivative, P * 1.2 X ICf 19 s s" 1 , 2,3 the pulsar is losing energy at a 
rate £ = 1.2 X IQ 36 erg s" 1 , only 2 X 10~ 3 times that of the Crab pulsar, 

Even if all of this energy loss is gravitational radiation, the star must be 
close to axi symmetry, with an eccentricity e = / 17 I45 < 4 X 10" 9 , 

where is the quadrupole moment perpendicular to a and is the moment of 
inertia in units of 10 45 g cm 2 . Because quadrupole moments can also result 
from anisotropic pressure and magnetic field stresses 18 , this is only an upper 
limit to rotationally induced axi symmetry . Regardless of the origin of the 
rapid rotation rate of the pulsar, it must presently have a < >Q cr -jt s ^ ,37 ^ 
irpG) 1 / 2 , from Eqn (1). Its mean density must therefore satisfy <p> > 2 X 10 14 
g cm“ 3 , in tne approximation of homogeneity. Some neutron star equations of 
state predict mean densities below this value for certain values of the 
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mass. Figure 1 shows models for neutron stars calculated from a variety of 
equations of state, ranging from very soft to very stiff dependences of 
pressure on density. The soft equations of state, including the pure neutron 
(R) and pion condensate cores require higher densities tc provide the 

pressure necessary to hold the star up against collapse. The stiff equations 
of state, however, allow lower densities which appear to be incompatible with 
the condition above for masses below 1.5 M 0 . For neutron stars 
with M < 0.5 Mg, which may not deform into tri axial ellipsoids, the limit on 
ft would not apply. If the millisecond pulsar has a mass in the range 0.5-1. 5 
Mg, then, the stiff equations of state can be ruled out. It would of course 
be interesting to investigate the equilibrium figures of rotating neutron 
stars with these equations of state to verify the above result. 

Several evolutionary scenarios have been proposed for this pulsar, and 
these models predict a whole class of short period pulsars with low magnetic 
fields. Alpar et. al. 5 suggest that neutron stars in binary systems may be 
spun up by Keplerian accretion disks. The minimum period is the period of the 
Keplerian orbit at the Alfven radius, so that 

P > 6 x ICf 4 Bg 6/7 R 6 18/7 (M/M Q r 5/7 fi- 7 3/7 s. , (2) 

where Bg is the surface dipole field in units of 10 8 G, Rg is the stellar 
radius in units of 10 6 cm., M is the mass and ft^y is the accretion rate in 
units of 10l 7 g s" 1 . Neutron stars with weak magnetic fields should 
therefore have the shortest periods. However, the star can only be spun up 
to the point of becoming non-axi symmetric, P eq = 1 X 10 -3 (M/Mq) - ^ 72 R 6 372 s. 
[using Eqn (1)], beyond which the excess kinetic energy of rotation will be 
released as gravitational radiation. This will be the case for neutron stars 
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with magnetic fields, 




V-j u t- *+■■>-* 
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< 1.8 Rg 5/4 (M/M 0 ) 1/4 fl,J /2 . 


(3) 


As the accretion continues, the star will remain at the equilibrium period 
P e q, emitting gravitational radiation at a rate equal to the rate of transfer 
of rotational kinetic energy from the disk, 


£ GRAV = Ifi eq^ 

= 4.4 x 1G 36 I 45 (M/M q } 1 / 2 R 6 ‘ 3/2 B 2/7 flj /7 er§ s 


(4) 


where I 45 is the moment of inertia, and the rate of spin-up $ is that given by 
Ghosh and Lamb^. 

Alternatively, the pulsar could have been born in isolation with a very 
weak magnetic field, as has been proposed by Brecher and Chanmugam 12 . In this 
case, if the initial spin were beyond the critical bifurcation frequency, the 
star would have quickly relaxed back to through gravitational 

radiation. 

Discovery of a shorter period pulsar could place further constraints on 
neutron star structure. For example, the two stiffest equations of state (MF 
and TI) could be eliminated completely for M > 0.5 M g by the existence of a 
pulsar with a period less than ~ 0.9 ms, and the Bethe-Johnson equation (BJ) 
would be restricted to either very low or very high masses. Since the BJ 
equation of state is currently believed to be the most realistic, the 
discovery of a pulsar with P < 1 ms. would have very interesting implications 
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for neutron star structure. The constraint on certain equations of state may 
have consequences for calculations, such as neutron star cooling times, which 
have been found to depend on the structure of the st<>r.^ 

I would like to thank 0. Kazanas and S. Tsuruta for discussion and 
comments which contributed to this work. 
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FIGURE CAPTION 


Figure 1 - Gravitational mass versus radius (from Baym and Pethick^) for s 
variety of neutron star equations of state. References are, beginning with 
the softest equations, Tr',ir (Ref. 15), R (Ref. 16), BJ (Ref. 17), MF (Ref. 18), 
T1 (Ref. 19). The hatched line denotes the mean density cutoff at <p> = 2 X 
10^ g cm“ 3 , imposed by the rotation of the millisecond pulsar. The dashed 
line? are density lower limits which would be set by shorter period pulsars. 




